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Abstract 

We investigate random density matrices obtained by partial tracing 
larger random pure states. We show that there is a strong connection 
between these random density matrices and the Wishart ensemble of ran- 
dom matrix theory. We provide asymptotic results on the behavior of 
the eigenvalues of random density matrices, including convergence of the 
empirical spectral measure. We also study the largest eigenvalue (almost 
sure convergence and fluctuations). 



1 Introduction 

Physicists and computer scientists working with finite size quantum systems are 
often interested in properties of typical states, such as entanglement, entropy, 
etc. In order to estimate such quantities, one has to endow the set of states 
(pure or mixed) with a certain probability measure and compute averages with 
respect to this measure. It has been known for a certain while now that there 
exists an "uniform" (in a way which will precised later) measure on the set 
£ n of pure states of size n. However, the situation is less simple when dealing 
with density matrices: there is no widely accepted candidate for a "canonical" 
measure on the set V n of all density matrices of size n. 

One may find in the literature two classes of probabilities on T> n : 

• the induced measures, where random density matrices are obtained by 
partial tracing a larger random pure state, 

• the metric measures, where the measure is the volume element associated 
to a particular distance on T> n . 

Depending on the physical nature of the model, one may choose different 
measures from one class or the other. In this work we study the measures of 
the first class. 

The induced measures were introduced by Braunstein [2] and studied later 
by Hall [5] , Zyczkowski and Sommers [T31 [13] . In the first part of this work 
we provide a rigorous construction of these measures. In the second part, we 
give some new explicit and recurrence formulas for the moments and we study 
the asymptotic behavior of the spectrum of such random density matrices. Our 
approach is based on the connection with the well-known theory of Wishart 
random matrices. 
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Our paper is organized as follows. In section [3 we recall the construction 
of the induced measures, adding mathematical rigor to the existing literature. 
Section [3] is devoted to recalling some results on the Wishart ensemble and 
making explicit the connection with random density matrices. We deduce the 
distribution of the eigenvalues and we study the moments. In Section[4]we study 
two models of large random density matrices, providing results on the behavior 
of the spectrum. A discussion of the results as well as ideas for generalizing our 
work are presented at the end of the paper. 

2 From pure states to density matrices 

We start by introducing and motivating the model of random density matrices 
that we consider. 

As explained in the Introduction, one would like to endow the set of density 
matrices on a complex Hilbert space Ti with a "natural" probability measure. 
It turns out that there is no straightforward way of doing this, so one has to 
make some additional hypothesis in order for a probability measure to stand 
out as the most natural one. Our approach here is based on the definition 
of a density matrix as it is usually understood in the theory of open quantum 
systems. We consider that the system described by the density matrix is coupled 
to an environment and that the compound system is in a random pure state. 
More precisely, we shall make two assumptions: 

(Al) The system Ti is coupled to an environment IC and the compound system 
Ti ® IC is in a pure state 

(A2) The pure state \tp) is uniformly distributed on the set of all pure states 
on Ti®IC. 

The first assumption is motivated by a large class of models considered in 
physics or quantum information theory. The general framework is provided by 
a system Ti (gi IC in a pure state, isolated from its environment. Suppose that 
one has access only to the sub-system Ti. This may happen for several different 
reasons: IC may be not accessible (e.g. Ti and K, are in distant galaxies) or IC 
may be too complicated to study (a heat bath or a noisy channel, for example). 
In these situations, it is natural to make the assumption (Al). Let us turn now 
to the second assumption. If one has no a priori information on the systems Ti 
and IC, it makes sense to suppose (A2). Moreover, it turns out that there exists 
an unique uniform probability measure on the set of pure states of given size, 
so we shall consider uniform pure states on the compound system. 

However, there are situations when one of the two hypotheses (Al) or (A2) 
is not physically motivated. For instance, when one has no knowledge of an en- 
vironment coupled to the system }C, there is no reason to suppose (Al). Instead, 
one should use other probability measures, such as the Bures measure (see the 
discussion in [13]). On the other hand, even if (Al) corresponds to the physical 
reality, one may have extra information on the system Ti or IC (or both). For 
example, it may be that the state of the environment IC has a particular form; 
thus, it makes no sense to assume (A2) and our model would not be adapted to 
such situations. 

In the next section, motivated by the assumption (A2), we shall construct 
the uniform measure on the set of pure states of given size. Then, by partial 
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tracing, we shall provide the probability which verifies the assumptions (Al) 
and (A2). 

2.1 The canonical probability measure on the pure states 

In quantum mechanics, a pure state is described by a norm one vector in a n- 
dimensional complex vector space TL. The phase of pure states is not determined, 
i.e. 

|e*V) = V(9eR (1) 

In order to make this definition rigorous, we introduce the following equiva- 
lence relation on 7i \ {0}: 

i~i/<S'3A€C* such that x = Ay. (2) 

Definition 1. A pure state \ip) is an element of the quotient space (7i\{0})/ ~. 
We denote by £ n the set of pure states of size n. 

As all complex Hilbert spaces are isomorphic to C™, the set £ n is the set of 
rays in C™. We endow £ n with the associated quotient topology and the Borel 
(T-ficld. We now turn to the construction of the uniform probability measure on 

As stated in the assumption (A2), the probability on £ n should be the most 
uniform one, as there is no a priori information on the state \tp). In particular, 
as there is no preferred basis of 7i, the uniform measure should be invariant by 
changes of bases. In our framework (Tt is a complex Hilbert space), changes of 
bases are provided by unitary applications. As a consequence, we ask that the 
uniform probability measure should be unitarily invariant. 

Definition 2. A measure v n on £ n is unitarily invariant if 

v n {UA) = v n {A), 
for all unitary U £ U(n) and for all Borel subset A C £ n . 

It turns out that this condition is strong enough to completely specify a 
measure on £ n , i.e. there is an unique unitarily invariant probability measure 
on £ n . This follows from a well-known result in probability theory regarding 
group actions (see [7]). Let us recall it here. 

Let G be a topological group acting on a topological space X. We call its 
action transitive if for all x, y G X, there exists g G G such that y = g ■ x and 
proper if for all g € G, the application X 3 x ^ g-x is proper, i.e. the pre- image 
of a compact set is compact. We then have the following 

Theorem 1 ([7 ). Let G be a topological group which acts transitively and 
properly on a topological space X . Suppose that both G and X are locally compact 
and separable. Then there exists an unique (up to a constant) measure v on X 
which is G-invariant. 

In order apply this result to our situation, we consider the action of the 
unitary group U{n) on the set £ n by left multiplication. We obtain the following 
proposition. 
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Proposition 1. The action ofU(n) on £ n is transitive and proper and thus 
there exists an unique unitarily invariant probability measure v n on £ n . 

Proof. First of all, notice that this action is well defined: the class \Uifj) does 
not depend on ip, but only on the class \ip); we say that the multiplication by 
an unitary is a class application. In order to show that the action is transitive, 
consider two classes and \tp} and an unitary U S U(n) such that Ui/j = (p 
(such an unitary always exists). It follows then that U\ip) — \<p). Finally, the 
action is compact, as the set £ n is compact and the multiplication applications 
are continuous. Thus, the action verifies the hypothesis of Theorem [TJ and 
as a consequence there is an unique unitarily invariant measure on £ n . More- 
over, given the compacity of £ n we can choose the measure of unit mass, which 
concludes the proof of the Proposition. □ 

Existence and unicity being settled, one would like to dispose of more con- 
crete descriptions on the distribution v n . It turns out that there are two ways 
of doing that. 

First of all, let us recall the definition of a complex Gaussian random variable. 
Let X and Y be two independent real Gaussian random variables of mean and 
variance 1/2. Then Z = X + iY is said to have a complex Gaussian distribution 
of mean and variance 1. We denote by A/"c(0, 1) the law of Z. A complex 
vector (Z\, . . . , Z n ) is said to have distribution Af£(0, I„) if the random variables 
Z\, . . . , Z n are independent and have distribution A/c(0, 1). 

Consider now a complex Gaussian vector X ~ A/J?(0, I n ) and the projection 
application 



It is well-known in probability theory that the law of X is unitarily invariant in 
C™. This property remains valid for the projection H(X) and thus the law of \X) 
is unitarily invariant on £ n . As there is only one unitarily invariant distribution 
on £ n , we have \X) ~ y n . 

We can also obtain the law v n from another well-known probability distribu- 
tion, the Haar measure onW(n). In order to do this, consider a Haar-distributed 
unitary matrix U. Obviously, the distribution of U is unitarily invariant; the 
same will hold true for the first column Y of U and for its class \Y). Thus \Y) 
has distribution v n . We sum up these results in the following 

Proposition 2. 1. Let X be a random complex vector of law J\f^(0,I n ). 
Then the class \X) of X is distributed along v n . 

2. Let U be a random unitary matrix distributed along the Haar measure 
on U(n) and let Y be the first column of U. Then the class \Y) has 
distribution v n . 

2.2 The induced measure on density matrices 

In this section we effectively construct the induced measures on density matrices 
that will be studied in the rest of the article. As stated in the Introduction, the 
induced measure of parameters n and k is obtained as follows: 



II : C 



H — > £„ 

X i— ► \x) 



(3) 
(4) 
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• Consider a product space Ti ® fC of two complex Hilbert spaces H (of 
dimension n) and JC - the environment - of dimension k. This is the global 
space system + environment. 

• Take an uniform random pure state \ip) on Ji ® JC (see the assumption 
(A2)). 

• Consider the (random) pure density matrix \ip)(ip\ corresponding to the 
pure state \ip). 

• Take p — Trfc(\ip)(ip\), the partial trace of \ip){ip\ with respect to the 
environment IC. The law of the random variable p is the desired probability 
measure, which we shall note p n .k- 

As in our formalism \ip) is an equivalence class, we shall define the pure 
density matrix \ip)(ip\ by: 

ww = iv|^) e ^ (c) - (5) 

Clearly, ip i— > IV'X'^I is a class function (it does not depend on the repre- 
sentant ip chosen, but only on the class IV'))) so IV'XV'I is well-defined. The 
normalizing factor Tr(ip • ip*) appears because we want the matrix |'0)(V , | to be 
trace one; this could have been avoided by considering a norm one vector ip, 
since Tr(^ • ip*) = \\ipf- 

We now turn to the third step of the above construction and recall that the 
partial trace is the unique application Tr^c : M n k(C) — > M n (C) such that 

Tr((A(g>I K )B) = Tr(ATT K (B)) VAe M„(C),B e M n k{C). (6) 

Its expression for elementary matrices (ai, a 2 € H, bi, b 2 € JC) is 

Trx:[(ai ® hi) • (a 2 ® b 2 )*] = (b 2 , h) ■ a ia * 2 . (7) 

We have now all the elements needed for the definition of the induced mea- 
sures: 

Definition 3. The induced measure of parameters n and k is the distribution 
Pn.k of the random density matrix 

p = Tr K (|^)(^|), (8) 

where \ip) is an uniform pure state on 7i £g> IC of distribution v n k- 

In order to get a better understanding of the measure p n ,k, we write ip in 
an orthonormal basis {e^ ® f 3 ■; 1 < i < n, 1 < j < k} of H ® JC: 

n k 

*p = j2J2^ ei ®h- ( 9 ) 
i=i j=i 

Thus the matrix \ip)(tp\ has coordinates (in the same basis): 



tpij 4>i>j> 

Eo=lE^=l I^M 



www = ^r j ., l2 - (io) 
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After taking the partial trace, we obtain 



P«' - ™ — — r, — [2- t 11 ) 

Now, if we arrange the coordinates ipij of ip in a n x fc matrix X such that 
X(i,j) = ipij, we have 

X-A* 

P= — ; r< (12) 

' Tr(A -X*) y ' 

Several important remarks should be made at this point. First of all, consider 
U G U(n) and the density matrix p' obtained by replacing ip with (U ® Iac)"0 : 

p' = Tr K (\(U ®I k )tP)((U ®I,c)tP\). (13) 

By the properties of the partial trace, we have that p 1 = UpU*. But recall that 
the law of \ip) is unitarily invariant; it is thus invariant by U ® l/c (which is 
an element of U(nk)). Hence the law is invariant by unitary conjugation. 
Being positive, and thus self-adjoint, density matrices diagonalizc: 

p = VDV*, (14) 

with V an unitary and D a diagonal matrix with positive entries. The unitary 
invariance of corresponds to the fact that V is distributed along the Haar 
measure on U(n). Remains, of course, the question of the distribution of D, 
the diagonal matrix of eigenvalues, which will be answered in Section 13.21 (see 
Proposition 0J. 

Another important question concerns the law of the matrix X. Recall that 
the coordinates of X are those of ip, rearranged in a n x k matrix. Since the 
pure state \ip) is distributed along the uniform measure v n k, we know, by the 
second point of Proposition[2j that we can take for ip a complex Gaussian vector 
in C nk . Thus, the elements of X are independent, complex Gaussian random 
variables. 

Lemma 1. Let X be a n x k complex matrix such that the entries are inde- 
pendent identically distributed (i.i.d.) A/"c(0, 1) random variables. Then, the 
matrix 

X-X* 

P ~ Tt{X ■ X*) ^ ' 

has distribution p n ^. 

This lemma motivates the study of matrices of type W — X ■ X* , which will 
be taken up in the next section. 



3 Wishart matrices. Results at fixed size 
3.1 The Wishart ensemble 

This section is devoted to introducing the Wishart ensemble of random matri- 
ces. Introduced in the 1930's to study covariance matrices in statistics, Wishart 
matrices have found many applications, both theoretical (random matrix the- 
ory) and practical: principal component analysis, engineering, etc. Let us start 
by recalling the definition of the Wishart ensemble: 
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Definition 4. Let X be a n x k complex matrix such that the entries are i.i.d. 
A/"c(0, 1) random variables. The n x n matrix W — X ■ X* is called a Wishart 
random matrix of parameters n and k. 

In virtue of Lemma [TJ there is a strong connection between the distribution 
of Wishart matrices and the random density matrices we study. More precisely, 
if W is a Wishart matrix, then 

W 

(16) 



TrW 



has distribution ji nt k- 

We shall give a list of results on Wishart matrices that will be used later 
in the study of random density matrices. As the results are rather classical in 
random matrix theory, we will not supply proofs, but only references to the 
original papers. 

We start with a result on the eigenvalues of a Wishart matrix. Being of the 
form W — X ■ X*, Wishart matrices are positive and thus they admit n non- 
negative eigenvalues Ai, . . . , A n . The next proposition provides the distribution 
of the random vector (Ai, . . . , A„) (see [5]). 

Proposition 3. Let W be a random n x n Wishart matrix with parameters n 
and k. Then the distribution of the eigenvalues (Ai, . . . , A„) has a density with 
respect to the Lebesgue measure on R™ which is given by 

/ n \ n 

$W k (Ai , .... An) = CZ exp - J2 ^ II A- - A(A) 2 , (17) 



where C^ k is the constant 



n-= r(n+i-i)r(fe-j) 



-i 



and 



A(A)= J] (A, -A,). (18) 

l<i<j<n 

When studying large random matrices, one important question is to what 
resembles the spectrum of a random matrix in the limit n — > oo? In order to 
answer such a question, one introduces the empirical spectral measure 

1 " 

L n (W) = -y^Sx,, (19) 

i=l 

which is a random probability measure (it depends on W, which is random). 
It turns out that, almost surely, the random measures L n (W) converge to a 
deterministic probability measure, the Marchenko-Pastur distribution. 

Definition 5. For c s]0, oo[, we denote by /i c the Marchenko-Pastur probability 
measure given by the equation 



[i c = max{l - c, 0}S H — l[ a>b ](x)dx, (20) 



where a = (a/c — l) 2 and b = {*Jc + l) 2 
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The result is contained in the following theorem (see [1]). 



Theorem 2. Assume that c £E]0, oo[ 7 and let (k(n)) n be a sequence of integers 
such that lim^^oo k(n)/n — c. Consider a sequence of random matrices (W n ) n 
such that for all n, W n is a Wishart matrix of parameters n and k(n). Define 
the renormalized empirical eigenvalue distribution of W n by 

1 " 

i=l 

where Xi(W n ), • • • , X n (W n ) are the eigenvalues ofW n . Then, almost surely, the 
sequence (L n ) n converges weakly to the Marchenko-Pastur distribution p c . 

Another object of interest in random matrix theory is the largest eigenvalue 
of a large random matrix. The next result shows that in the Wishart case, 
it converges almost surely to the right edge of the support of the Marchenko- 
Pastur distribution; similarly to the Central Limit Theorem, the nature of the 
fluctuations is known (see [T] and [H]). 

Theorem 3. Assume that c £E]0, oo[ 7 and let {k{n)) n be a sequence of integers 
such that linin^oo k{n)/n — c. Consider a sequence of random matrices (W n ) n 
such that for all n, W n is a Wishart matrix of parameters n and k(n), and let 
X m ax{W n ) be the largest eigenvalue ofW n . Then, almost surely, 

lim -X max (W n ) = (V^+l) 2 - (21) 

n — >od n 

Moreover, the following limit holds in distribution 

n^oo nV3(l + a/c)(1 + 1/y/c) 1 / 3 

Here, W2 is the Tracy- Widom law of order 2; as even the definition of this 
probability distribution is well beyond the scope of this work, we encourage the 
reader to look it up in [14] , the original paper of Tracy and Widom. 



3.2 The spectrum of a density matrix 

Recall from Section 12.21 that when considering the diagonalization of a random 
density matrix 

p=VDV*, (23) 

the unitary matrix V is distributed along the Haar measure on the unitary group 
lA(n). In this section we compute the distribution of the diagonal matrix D, i.e. 
the spectrum of a density matrix with distribution fi„ t k- 

Here, as well as in the next section, the parameters n and k will be fixed, 
and we shall assume that k > n. If n > k, by a property of the partial trace 
application, the matrix p will have n — k null eigenvalues and k eigenvalues 
identical to those of the density matrix 

o- = Tr n (\i>)(<p\). 

In consequence, the study of the spectrum of p is equivalent to the study of 
the spectrum of a. Moreover, the size of cr's environment (n) is larger than 



8 



the dimension of a itself (k), and we can apply the first case. In conclusion, 
whenever n is larger than k, we interchange n and k, and we append n — k null 
eigenvalues to the spectrum of p. 

Recall that if W is a Wishart matrix of parameters n and k, then p = 
Wj Tr W has distribution fi n ,k- It follows that if (Ai, . . . , A„) are the eigenvalues 
of W and (Ai, . . . , A„) are those of p, then we have 

A, = , , Vl<i<n. (24) 

2^=1 Aj 

As the trace of a density matrix equals one, the (random) vector (Ai, . . . , A„) 
is confined in the (?i— l)-dimensional probability simplex — {(xi, ■ ■ ■ ,x n ) G 
MIL : Yli=i x i ~ •"■}■ Note that A„ is a function of Ai, . . . , X n -i, so we will show 
that the distribution of (Ai, . . . , A n _i) admits a density w.r.t. the Lebesgue 
measure on X n _i. 

Proposition 4. The distribution of the (unordered) eigenvalues Ai(p), . . . , A„_i(p) 
has a density with respect to the Lebesgue measure on £„_i given by 

n 

* w ,fc(Ai, • , A„-i) = C n , fe n (A,) fe -"A(A) 2 , (25) 

i=l 

n r(nfc) 

n j= o r(n + l-j)r(fc-j) 



Remark 1. In i/ie formula h25}) , there are only n — 1 variables; X n is noi a 
variable, but merely the notation X n = 1 — (Ai + ■ ■ • + A n _i). 

Proof. Let us start from the Wishart distribution of eigenvalues and consider 
the change of variables 

(Ai, . . . , A„) h-> (Ai, . . . , A„_i, S) (27) 
h- (Ai/5, . . . , A n _i/5, 5) = (A x , . . . , An_i, S). (28) 

where 5 = X)"=i Aj is the sum of the Wishart eigenvalues. The Jacobian of this 
transformation equals 1/S n , and we get 

n 

$ { ~ X f(\i, . . . , A,,-!, S) = C™ k exp(-S) Y[ (SXi) k ~ n A(SX) 2 S n ^ 1 . (29) 

i=l 

We get now to the crucial point of the proof. Clearly, the above density factorizes 
as 



$5 S) (A 1 ,...,A n _ 1) 5) = C> v fe x 



nA.f-A(A) 2 



i=l 



x [5 ,nfe - 1 exp(-5)] . (30) 



Hence, the normalized eigenvalues (Ai, . . . , A„_x) and the sum of the Wishart 
eigenvalues S are independent random variables. 
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In order to compute the distribution of (Ai, . . . , A„_i), it suffices to take the 
marginal with respect to S; we get 



*„ )fc (Ai, . . . , A„_i) = C„, fc J] Af-"A(A) 2 , 



(31) 



where 



poo 

,,,.=< „.,,.- / S nk - l e- s dS = T(nk)C™ k = 
Jo 

T(nk) 



uUnn+i-mk-j) 



(32) 
(33) 
□ 



As a byproduct of the proof, we also obtain the following characterisation 
of the induced measure. 

Corollary 1. The law of a random density matrix p of parameters n and k is 
the law of a Wishart matrix W of the same parameters conditioned by Tr W = 1 . 

Proof. From the formula ([3H|) we see that the normalized eigenvalues and the 
trace of a Wishart matrix are independent random variables. Thus, taking the 
marginal with respect to the trace is equivalent to conditioning on the event 
(TrVT = 1). Note however that (Tr W = 1) has zero probability. □ 

In the Figure [T] we have plotted the density functions for n = 2 and several 
values of k using the analytic formula (f25|) . For n — 3 we have randomly 
generated random density matrices and plotted the probability simplex E2 along 
with the points corresponding to the spectra (Figure H|) . We notice that for 
large values of k (the size of the environment), the spectrum concentrates to the 
middle point in S„_i. This is a general phenomenon and it will be studied in 
section 14.11 





Figure 1: Theoretical eigenvalue distribution for (n = 2, k = 2), [n = 2, k = 10) 
and ( n = 2,k = 50). 



3.3 Moments 

The aim of this section is to provide formulas for the moments of order q of 
a random density matrix of distribution pL n ^. In order to do that, we shall 
introduce the some notation: E n) fc[-] will denote the expectation with respect 
to the law pi n ^ and the expectation with respect to the law of Wishart 
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Figure 2: Empirical eigenvalue distribution (5000 matrices) for (n = 3, k = 3), 
(n = 3, jfc = 10) and (n = 3, jfe = 100). 

matrices with parameters n and k. We will use the corresponding result on the 
Wishart ensemble and derive explicit formulas, as well as recurrence relations. 
The following proposition provides a bridge between the moments of a density 
matrix and those of a Wishart matrix with the same parameters. 

Proposition 5. Let E„ ! fc[Tr(p' 3 )] be the moment of a random density matrix of 
parameters n and k and let K^ k [Tr(W q )} be the moment of a Wishart matrix 
having the same parameters. Then, 

[Tr(W^)] 
nk(nk + 1) ■ ■ ■ (nk + q — 1) 
Proof. By using the same technique as in the proof of the Proposition H] we get 

< fe [Tr(W<0] = E», fc [Tr(^)] r( ^+ g) ; (35) 

which is the same as equation □ 

We can find in the literature different explicit and recurrence formulas for 
E™ k [Tr(W q )]. From the one in 0], we get 

M ^ ,J -r(n* + ? )^ ( 1} (?-i)!(j-l)! ' (36) 

where [a] q = a(a — 1) • • • (a — q + 1). The recurrence formula (see [1]) 

(,-»((, -!)»-(*-»)») Etiil[Tr(p ,- 2)] (38) 



(nk + q - l)(nk + q- 2)(q + 1) 
allows us to easily compute some averages: 

E„, fc [Tr(p 2 )]^-^, (39) 

rm , ■„ n 2 + 3nfc + k 2 + 1 . . 

E ^ Tr ^= (^TTPT2) ' (40) 

rm / 4m n 3 + 6n 2 fc + 6nfc 2 + £r 3 + 5n + 5/c 

E ".*F & » 4 )]= (n fc + l)(n fc + 2)(n fc + 3) ' ^ (41) 

These formulas are consistent with the ones of [13] and ITS]. 
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4 Asymptotics 

The last part of this paper is devoted to the study of random density matrices 
corresponding to large systems. We shall consider two models, both motivated 
physically: 

1. In the first model, the size of the density matrix n is constant and the size 
of the environment k tends to infinity. Such a situation arises typically 
when one studies a small system (a qubit, a pair of qubits, etc.) coupled 
to a much larger environment. We show that in the limit k — > oo, den- 
sity matrices distributed along yL n .k converge to the maximally mixed (or 
chaotic) state I jn. 

2. In the second model, both n and k tend to infinity and k/n — > c > 
0. This model describes a large system coupled to a large environment 
with constant ratio of size (dim K, / dim TL sa c) . In this case we show 
that the spectral measures of density matrices of law \i n ^ converge to a 
deterministic measure known in random matrix theory as the Marchenko- 
Pastur distribution (see Definition [5]) . We also study the convergence and 
the fluctuations of the largest eigenvalue of random density matrices. 



4.1 The first model 

Consider the density function of ji ny k with n fixed and k — > oo: 



*» )fc (Ai, . . . , An-i) = C n , k n Af-"A(A) 5 



(42) 



As n is fixed, the Vandcrmonde factor A(A) is constant. The other factor, 
properly normalized in order to get a probability density, is the Dirichlet measure 
of parameter a = k — n + 1 : 



$n,fc(^i> • • • ) -Vi-i) — C n ,fc J\ X " 



(43) 



The next result is well-known in probability theory. We shall sketch its proof 
for the sake of completenss. 

Theorem 4. The Dirichlet measure converges weakly as a — > oo to the Dirac 
measure £(i/ n ,...,i/n) 

Proof. The idea behind the proof is to show that the variance of a Dirichlet- 
distributed random variable converges to as its parameter converges to infinity. 
Let X be such a random variable. X has a density with respect to the Lebesgue 
measure on the probability simplex given by: 



f{x\ , . . . , x n ) 



It is easy to compute 
E 











2 






X - 








= nE 



T(na) 
T(a) n 



2xi 
n 



n 



a + 1 



1 



na + 1 n 



0. (44) 
□ 
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As the maximally mixed state I /n is the unique state having spectrum 
{1/n, 1/n}, we get: 

Corollary 2. Density matrices of the first model converge almost surely to the 
maximally mixed (or chaotic) state I /n. 

Remark 2. The same result can be obtained by an entropic argument. It turns 
out that the mean von Neumann entropy S(p) — — Tr(/9log/o) can be computed 
for a random density matrix distributed along \x n ^: 

EmF(p)]= E^- 

i=k+l 

This formula has been conjectured by Page and has been subsequently proved 
(see using various methods. Let us explain how it implies the corollary. 

First, fix n and let k grow to infinity, as in our model. The mean entropy is 
easily seen to converge to logn. This turns out to be the maximum von Neumann 
entropy for a system with n degrees of freedom. It is attained at the state I/n, 
the unique state of maximum uncertainty. 



4.2 The second model 

In the second model, both the size of the density matrix and the size of the 
environment tend to infinity. In order to use the results on the Wishart ensemble 
(Theorems [5] and [3]), we need appropriate results on the behavior of the trace S 
of a Wishart matrix. 

Lemma 2. Assume that c G]0, oo[, and let (k(n)) n be a sequence of integers 
such that lmin^oo k(n)/n = c. Consider a sequence of random matrices (W n ) n 
such that for all n, W n is a Wishart matrix of parameters n and k(n). Let 
S n = Tr W n be the trace of W n . Then 

S 

— — — — > 1 almost surely (45) 
nk(n) 

and 

Sn ; nk{n Uu(0M (46) 
yj nk(n) 

where " =>■ " denotes the convergence in distribution. 

Proof. Recall that W n = X n ■ X*, when X n is a n x k(n) matrix with i.i.d. 
complex Gaussian entries. We have 

n k{n) n k(n) 

s n = E E i*«i 2 = E E( Rc (^) 2 + MX,,) 2 ). (47) 
i=i j=i i=i j=i 

The random variables {Re(Xy), lm(Xij)}ij are i.i.d. with distribution Af(0, 1/2) 
and thus, by the law of large numbers, we have, almost surely, 

lim Trh = ( 48 ) 

n^oo 2nk(n) 2 

completing the proof of the first result. The second result follows from the 
Central Limit Theorem. □ 
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We can now state and prove the analogue of Theorem [5] for random density 
matrices: 

Theorem 5. Assume that c £]0, oo[, and let (fc(n)) n be a sequence of integers 
such that lim n _ ) . 00 k(n)/n — c. Consider a sequence of random density matrices 
(pn)n such that for all n, p n has distribution p n ,k(n)- Define the renormalized 
empirical distribution of p n by 



1 

L n = -^2S cnXi (p n ), (49) 

i=l 



where Xi(p n ), - • • ,X n (p n ) are the eigenvalues of p n . Then, almost surely, the 
sequence (L n ) n converges weakly to the Marchenko-Pastur distribution p c . 

Proof. We know (Theorem [5]) that the empirical distribution of eigenvalues for 
the Wishart ensemble 

n 

n * — ' 

i—l 

converges almost surely to the Marchenko-Pastur distribution of parameter c. 
Recall that the eigenvalues of the density matrix p n = W n / Tr(W n ) are those 
of W n divided by the trace S n of W n . We have thus the following formula for 
the empirical spectral measure of p: 

~ n 1 n 

Ln = ~^Z 5 cn\^W n )/S n = -XX-U,(W„).^- ( 51 ) 

i=l i=l 

The last equation is the same as equation (|50|) with the Dirac measures per- 
turbed by a factor of cn 2 / S n which converges, almost surely, to 1 (by the pre- 
ceding lemma). We are now going to show that such a perturbation does not 
change the limit in distribution. In order to achieve this, recall that when the 
limit measure is compactly supported, the convergence in distribution is equiv- 
alent to the convergence of moments. If we compute the g-th moment of the 
measures and L n , we find: 

(x\L™) = ~f2(n- 1 \ i (W n )) q , (52) 



n 



and, respectively, 



(Ai«)4E(^W) ? ' f • (53) 

i—l ^ ' 

These expressions have the same limit as n — > oo for all q, and thus L n converges 
to the Marchenko-Pastur distribution. □ 

In the Figure El we have plotted for several values of c and large n and k a 
histogram of the spectrum for one density matrix and the theoretical density 
of the Marchenko-Pastur distribution (see Remark [3]). We can see that the 
empirical histogram matches closely the theoretical curve for rather mild values 
of n (here n = 1000). 

We now turn to the study of the largest eigenvalue of random density ma- 
trices. As before, we use the corresponding result on the Wishart ensemble 
(Theorem [3]) and the control over the trace (Lemma [2]): 
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Figure 3: Empirical and limit measures for (n = 1000, k = 1000), (n = 1000, k = 
2000) and (n = 1000, k = 10000). 



Theorem 6. Assume that c s]0,oo[, and let (k(n)) n be a sequence of integers 
such that linifj^oo k(n)/n = c. Consider a sequence of random matrices {p n )n 
such that for all n, p n has distribution /x n ,fc( n ), and let X max (p n ) be the largest 
eigenvalue of p n . Then, almost surely, 

lim cnX max (p n ) = (\/c+ l) 2 . (54) 

n^oo 

Moreover, 

n 2 / 3 \cnX max (p n ) - (Vc+ 1) 2 1 
lim t =— , _ : ,„ = W2 in distribution. (55) 

(1 + V^)(i + i/\/H) 1/3 



n — >oo 



Proof. By the first part of Theorem [3J the (normalized) largest eigenvalue 
h^max(Wn) of a Wishart matrix converges almost surely to (yfc + l) 2 . Ob- 
viously, we have 

( W n) fK ^ 

■A-max(Pn) — 5 , \00 ) 

"11 

and, by the Lemma S n /(cn 2 ) converges (almost surely) to 1. Finally, we 
obtain formula (j5"4|) . 

For the second part of the theorem, what we need to do, normalizations 
apart, is to show that the trace of a Wishart matrix fluctuates less than the 
largest eigenvalue. For the Wishart case, we have 

Xmax(W n ) = n(V~c + l) 2 + n 1/3 (l + V~c)(l + l/V~c) 1/3 (m + o(l)), (57) 

and 

S n = nk{n) + ^ftMn){N + o{l)). (58) 

Again, X max (p n ) = X max (W n )/S n and after simplifications, one obtains the 
desired formula ([55"|) . □ 

Remark 3. Note that Theorem^ and the first part of Theorem^ deal with al- 
most sure convergences. This means that when considering sequences of random 
density matrices of increasing size, the respective convergences will fail only on 
a set of null measure. This is to be compared with typicality results for random 
density matrices obtained recently in fl(J\/ by concentration of measure tech- 
niques. These results give bounds (at fixed matrix size) on the probability that 
a random matrix is far from its expected value, while our results deal with the 
more subtle convergence of rescaled quantities, such as the spectral distribution 
or the largest eigenvalue. 
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5 Conclusions 



We investigated random density matrices distributed along the so-called induced 
measures. After introducing them as partial traces of larger random pure states, 
we provided some explicit and recurrence relations for the moments of such 
density matrices. Using results on Wishart matrices, we then considered large 
density matrices. In a first model, a fixed size system was coupled to a very 
large environment; we showed that an uniform pure state on the compound 
system corresponds to the maximally mixed (or chaotic) density matrix on the 
fixed-size system. In parallel with Wishart matrices, we studied the regime 
dim/C/ dim7i — + c. We obtained the almost sure convergence of the empirical 
spectral measure and of the largest eigenvalue, as well as the fluctuations of 
the largest eigenvalue. Results from random matrix theory were easily adapted 
for density matrices. Other important quantities, such as correlation functions, 
require a more detailed analysis, and this will be the subject of further work. 
Also, it may be interesting to study such asymptotics for other probability 
measures on density matrices, such as the Bures measure. 

Acknowledgment: The author would like to thank Guillaume Aubrun for 
useful ideas which led to several simplifications in some proofs. 
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